Output from an atom laser: theory vs. experiment 
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Atom lasers based on rf-outcoupling can be described by a 
set of coupled generalized Gross-Pitaevskii equations (GPE) . 
We compare the theoretical predictions obtained by numeri- 
cally integrating the time-dependent GPE of an effective one- 
dimensional model with recently measured experimental data 
for the F=2 and F=l states of Rb-87. We conclude that the 
output of a rf-atom laser can be well described by this model. 



The discovery of Bose-Einstein condensation (BEC) in 
trapped atomic gases |l|^| and the subsequent demon- 
stration of macroscopic coherence of Bose condensed 
gases HH have raised a lot of interest in creating coher- 
ent atomic beams using such condensates. Experimental 
setups based on Bose-Einstein condensates that produce 
coherent atomic beams are now called atom lasers. The 
first experimental realization || was based on a pulsed 
outcoupler; strong radio- frequency pulses were used to 
flip magnetically trapped atoms to untrapped states such 
that they could leak out of the trap. The method of rf- 
outcoupling has been extensively investigated theoreti- 
cally ||7-|L3|) whereas the possibility of using Raman tran- 
sitions |Q has not been treated in such detail. Most 
recently, the outcoupling process was treated taking also 
thermal excitations within the trap into account |15| . 

Experiments have also made some remarkable prog- 
ress. In [pi a Bose condensate is released from a mag- 
netic trap into an optical standing wave leading to mode- 
locked like pulses of atoms falling downwards. A method 
using Raman pulses is demonstrated in [ p.7| . There, the 
outcoupled atoms are moving horizontally and are pulsed 
with a very high frequency resulting in a quasi-continuous 
beam. Finally, Bloch et al. [Q use a rf-outcoupler within 
a very stable magnetic setup, allowing to operate the 
atom laser in a regime close to the weak-coupling regime 
dealt with in [n0[. As in all other cases there is no refill 
mechanism for the Bose condensate in the trap that bal- 
ances the outcoupling loss. Nevertheless, the outcoupling 
rate is so low that the output is similar to a cw-laser with 
a slowly decreasing intensity. 

In the following, we want to model the situation in |L8| . 
The trap is elongated horizontally. As the atomic beam 
is directed downward, we approximate the trap with a 
ID model in the vertical z-direction taking account of 
the other two dimensions in a way explained in the next 
section. 



I. THEORY OF COUPLED GROSS-PITAEVSKII 
EQUATIONS 

Condensates of weakly interacting Bose gases at zero 
temperature are well described by the Gross-Pitaevskii 
equation (see |19| and references therein). As it is a 
nonlinear Schrodinger equation for a macroscopic wave- 
function it not only accounts for the density proper- 
ties but also for the coherence inherent in Bose con- 
densates. In a generalization of the GPE was in- 
troduced to treat a system of trapped and untrapped 
states of a Bose gas coupled via an rf field with a fre- 
quency w r f. Here, we want to model an atom laser for 
Rb-87 in its F — 2 or F = 1 hyperfine-manifold. The 
2F + 1 Zccman sublevels are represented by macroscopic 
wavefunctions ip m (m € {— F, . . . , F}) so that the system 
of generalized GPEs after applying the transformation 
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Thereby we have assumed that all Zeemann levels inter- 
act with the same s-wave scattering length ag = 110 asohr 
so that U = A-K% 2 a N/M. The total density devided 
by the number N of atoms in the system is denoted by 
||V>(r,£)|| 2 = J2 m \^m{r,t)\ 2 (we normalize to 1, not to 
N). The Rabi frequency is determined by the strength 
2? r f of the rf-field and given by hQ = gp^o^B-^/2. It 
is modified by the matrix elements of the angular oper- 
ators F± represented by c m = \/F{F + 1) — m(m + 1). 
For F = 1, this results in c m = v2; for F = 2, c m takes 
values 2 and \/6 depending on m so the Rabi frequencies 
differ for different m — > m! transitions. 

Ioffe-type traps often used in experiments are usually 
elongated in the horizontal plane being axial symetric in 
the remaining directions. As mentioned above, we want 
to model an atom laser based on such traps by ID GPEs 
where the axis under consideration is the vertical axis of 
the trap. We choose coordinates where the z-axis points 
downwards, the long axis of the trap is denoted by y, the 
short horizontal one by x. Taking gravity into account, 



the total effective potentials in z-direction then are given 
by 

V m , e s(z,t) = sgn(g F )mMu; 2 z 2 /2 + mhA 

-Mgz + g 1B Mz,t)\\ 2 . (1.2) 

A = 7iu; r f — V a g is simply the detuning of the transitions 
(VoB = -ffFMBohr-Bog, w r f > for F = 1 and Wrf < for 
F = 2). 

The effective interaction strength <?id is determined in 
such a way that the chemical potential in Thomas-Fermi 
approximation in the ID model equals that for the full 
3D situation for TV particles in a trapped state TO tra p ■ It 
turns out to be equal to 
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with a z = y / h/(Mui z ). The ID version of Eq. (1.1) now 
reads 

d ( h 2 d 2 \ 

iTi—il) m (z,t) = I -— —-2 + lWff(M)J i/} m (z,t) 

+fj£l};(c m >8 mjm ' + i +Cm5 m ,m'-l) 1pm' (z,t). (1.4) 



Gravitation leads to a sag of the trapped atomic cloud. 
The shift of the minimum of the potentials away from 
the trap center at z — depends on to and is given by 
Zm.sne = g/(\ m \u 2 ). The effective potentials always cross 

Usually, only the positive 
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= ±a zy /2A/(hu z ). 
sign plays a role because the negative value is outside 
the condensate (cf. inset in Fig. P)). 
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FIG. 1. Plot of the densities for m = -1,0 (F = 1) in 
the trap for weak outcoupling. The trapped state (thick line, 
scaled down by a factor 20) is only disturbed slightly, whereas 
the untrapped state (thin line) can be approximated by an 
Airy function (dotted line shows the magnitude). The inset 
shows the effective potentials. 



For weak coupling fl < lu z the trapped atoms are 
only disturbed slightly by the outcoupling process p0| , 
after a certain propagation time the wavefunctions are 
quasi-stationary. In Fig. [j] we plot the densities of the 
to = —1,0 states for F = 1 in such a situation. The 
outcoupling takes place primarily at the resonance point 
z Ics . The untrapped state can be approximated by an 
Airy function, its strength being determined by the out- 
coupling rate which is known analytically ptj| , p0[ . 

Why is the location of maximal outcoupling identical 
to z rcs ? This is not so obvious; if one considers the same 
situation without the mean-held contribution in V mte s, 
the density maximum of the untrapped state is located 
at the classical turning point in its potential, which de- 
pends on the energy (analogous to photo-dissociation of 
molecul es) . If the mean-field contribution is present like 
in Eq. (1.4) the classical turning point is shifted to the 
resonance point z res (see inset in Fig. (0)), and hence the 
point of maximum density coincides with z ros . 



II. NUMERICAL RESULTS 

In this section we present some numerical results ob- 
tained by propagating the set of coupled GPEs Eq. ( |l.4| ) 
and compare them to experimental data for the same set 
of parameters. The time propagation was accomplished 
by a usual split-operator method using FFT on a ID grid. 
The starting condition is always the condensate ground 
state obtained by imaginary time propagation. 

The basic parameters in our calculations are taken 
from ]lq| , the \F = 2, to = 2) -state is trapped with 
frequencies u>± = 2tt x 180 Hz and wm = 2tt x 19 Hz. 
Throughout the paper we use harmonic oszillator units 
with respect to the \m\ = 1 trapping frequency for 
the vertical motion Wtrap = 27r x 127 Hz which hap- 
pens to be the same for both F = 1 and F = 2 
(<7f=i = — 1/2, <7f=2 = 1/2). Thus, the length unit is 
a z = 0.95 /im, time is measured in 1/wtrap = 1.3 ms. 

The Rabi frequencies used in |18[ are all bigger than 
Wtrapi so we are not in a weak coupling regime. During 
the outcoupling process, the condensate is depleted sig- 
nificantly. In Fig. (|2|), one can see how fast this takes 
place. The populations of the states show some small 
oscillations that are reminiscent of the Rabi oscillations 
occuring at very strong RF fields where there is hardly 
no outcoupling [H,13|. Such oscillations can be seen in 
experiment as a sort of pulses pl[ . 
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FIG. 2. Time evolution of the populations for F = 1 atom 
laser. The condensate in the m — — 1-state is depleted very 
fast (fi/^trap ~ 3-5), the time axis extends over 12.5 ms. The 
period of the small oscillations corresponds to the Rabi fre- 
quency Q. 

In [lg] Bloch et. al. presented measurements of the 
number of particles remaining in the \F — 2, m = 2)- 
state after a fixed time of outcoupling depending both on 
the strength and the frequency of the radio freqency. Fig- 
ures (H) and (|4|) display our numerical data for their pa- 
rameters and the experimental values. In the calculations 
we took only the m = 2,1,0-states into account. Tests 
with more states show that this is enough, as the missing 
m = — 1, — 2-states are populated only very weakly. 
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FIG. 3. Number of atoms remaining in F — 2, m = 2 af- 
ter 20 ms of outcoupling. Initially there are N — 7 x 10° 
atoms. The detuning A is tuned to maximal outcoupling, so 
Zrcs ~ -22, sag- B r { is varied from 0.1 to 1.0 mG leading to Rabi 
frequencies from 0.44 to 4.4 kHz. Lines are drawn for guiding 
the eye only. 




FIG. 4. Same 

in Fig. (H), but now Q, = 2ty x 700 Hz = 5.5oj t ra P is kept 
fixed and the detuning A (resp. radio frequency) is varied. 
The experimental data is shifted in frequency to match the 
minimum of the theoretical curve; this is necessary because 
the offset field B a of the trap is not known with high enough 
precision. Again, lines are there for guiding the eye. 

The theoretical values fit the experimental data at least 
qualitatively. The differences are most probably due to 
dimensional effects: in the ID model there is only a 
resonance point at z res , whereas in 3D the outcoupling 
takes place on a resonance surface defined by the dis- 
tance r res = z ros to the trap center. Accordingly, the 
outcoupling rate in 3D depends on the density of the 
trapped state averaged over the resonance surface which 
is smaller than the density on the z-axis. Thus, the out- 
coupling yield in 3D is smaller than in ID because there 
is only the "on-axis" density. 

This effect is even stronger for the F = 1 results shown 
in Fig. (0); the experimental data are again from the 
Hansen group |21|. In this case, all three Zeeman sub- 
levels were propagated because they are all significantly 
populated, (see e.g. Fig. (||)). We have tried to quantify 
the aforementioned difference between the ID- and 3D- 
situation including gravity at least in the weak co upling 
limit in the spirit of the perturbational rates from |10 20 
but have not found a conclusive answer. 
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FIG. 5. Same as in Fig. (§) but for F = l,m = -1. and 
f2 = 27r x 312 Hz = 2.4utra P - The curve at weaker coupling 
(fi = l.lcJtrap) clearly exhibits an asymmetry reflecting the 
density distribution of the trapped state (see text). Mind also 
the bump around A = 120 entrap in the "theory" curve. 

A closer look on the spectral data in Figs. (||) and (||) 
reveals two further details: from a weak coupling analysis 
one expects the curves in these figures to reflect directly 
the density distribution of the trapped states transfered 
to the frequency domain via z — a z W2AJ(favtrap) ■ This 
leads to an asymmetry around the minimum that can be 
seen in the curve at weaker coupling shown in Fig. (pi). 

The second detail concerns the small bump visible in 
the "theory" -curve in Fig. (||) around A = 120w tr ap- It 
occurs only for intermediate couplings. We have no ex- 
planation for this feature. 



III. CONCLUSIONS AND OUTLOOK 

In this paper, we have described a numerical analysis 
for an RF-atom laser. The output yield of such a de- 
vice can be determined in a satisfactory way though the 
restriction to one dimension does not allow for a exact 
quantitative comparison with the real world data form 
experiment. 

What is lacking is a more analytical treatment of the 
outcoupling process for intermediate coupling strengths, 
which would probably also explain the bump in Fig. (|^). 
Furthermore, it would be nice to have a at least 2D pic- 
ture for the outcoupled wavefunction to discuss things 
like the beam divergence. A 2D treatment would also be 
able to account for the long axis of the trap where collec- 
tive excitations of the condensate can occur much easier 
than in the other two directions (cf. p5|). 

Another intersting subject is an atom laser driven by 
two radio frequencies. This may lead to a pulsed oper- 
ation mode like in a mode-locked laser. We are investi- 
gating such setups, results will be published elsewhere. 
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